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(N ■ ABSTRACT 

bo: 

Radial-velocity planet search campaigns are now beginning to detect low-mass 
"Super-Earth" planets, with minimum masses M sin i< 10M®. Using two 
', independently-developed methods, we have derived detection limits from nearly 

four years of the highest-precision data on 24 bright, stable stars from the Anglo- 
Oh \ Australian Planet Search. Both methods are more conservative than a human 

analysing an individual observed data set, as is demonstrated by the fact that 
both techniques would detect the radial velocity signals announced as exoplanets 
for the 61 Vir system in 50% of trials. There are modest differences between the 
methods which can be recognised as arising from particular criteria that they 
adopt. What both processes deliver is a quantitative selection process such that 
one can use them to draw quantitative conclusions about planetary frequency 
and orbital parameter distribution from a given data set. Averaging over all 24 
CN \ stars, in the period range P <300 days and the eccentricity range 0.0 < e < 0.6, 

^ \ we could have detected 99% of planets with velocity amplitudes K>7.1 ms _1 . 

\ For the best stars in the sample, we are able to detect or exclude planets with 

K>3 ms _1 , corresponding to minimum masses of 8M® (P=5 days) or 17M ffi 
O \ (P=50 days). Our results indicate that the observed "period valley," a lack of 

giant planets (M > 100 M®) with periods between 10-100 days, is indeed real. 
| However, for planets in the mass range 10-100 M e , our results suggest that the 

^ ' deficit of such planets may be a result of selection effects. 

Subject headings: planetary systems - techniques: radial velocities 



O 



department of Astrophysics, School of Physics, University of NSW, 2052, Australia 
2 Anglo- Australian Observatory, PO Box 296, Epping, 1710, Australia 

3 Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, Herts ALIO 9AB, 
UK 

4 Department of Terrestrial Magnetism, Carnegie Institution of Washington, 5241 Broad Branch Road, 
NW, Washington, DC 20015-1305, USA 

5 Faculty of Sciences, University of Southern Queensland, Toowoomba, Queensland 4350, Australia 



-2- 



Introduction 



Extrasolar planet detections are inexorably pushing toward ever-lower masses, and it 
seems that each new mass threshold spawns an unavoidably anthropocentric moniker. First, 
in 1995, were the "hot J up iters" (or "ro asters"; Sudarsky et al . 2000), with minimum masses 
m sin i greater than the mass of Jupiter ([Mayor fc Quelozlll995h . The earliest planet discover- 



ies fell into this category due to the well-known radial- velocity selection bias toward detecting 
massive planets in short-perio d orbits (P < 10 days) . The discovery of several Neptune- mass 



17.4 M e ) planets in 2004 (IMcArthur et all 12004 : ISantos et all 12004 ; Bonfils et all 120051 ) 



brought the "hot Neptune" into common parlance. Most recently, the discovery of planets 
with m sin % less than 10 Earth masses has prompted the somewhat misleading appellation 
of "super-Earth." These planets, of course, are not likely to bear any resemblance to Earth, 



Adams et all 120081 ) 



almost all of 


these 


(Sealer et al. 


2007; 



The Anglo-Australian Planet Search (AAPS) has been in operation at the 3.9m Anglo- 
Australian Telescope for 11 years, and has published 36 exoplanet discoveries. By integrating 
for at least 20 minute s on each target in order to average over the stellar p-mode oscillations 
(jO'Toole et al.l 120081 ) , the AAPS now ach ieves velocity precision approaching 1 ms _1 for 
bright, inactive stars (jO'Toole et al.ll2009d ). Since 2005 June, the AAPS has been routinely 
implementing this observing strategy, which has (together with long continuous observing 
blocks) proven instrume ntal in our recent discovery of a super-Earth and two Neptune-mass 
planets orbiting 61 Vir (jVogt et al.l 120101 ) . 



In lO 'Toole et al.l (j2009d ). we explored the frequency of planets with periods less than 
~20 days. In that work, we derived detection sensitivities from a continuous 48-night ob- 
serving campaign targeting 24 bright, stable stars. The stars have spectral types between 
GO and K5, and they r ange i n mass from 0.77-1.28 M Q ; complete details are given in Ta- 



ble 1 of lO 'Toole et al.l (l2009d ). In this paper, we consider periods ranging up to 300 days, 
now making use of nearly 4 years of high-quality data on these 24 stars from the AAPS. 
The observational data are summarised in Table 1, and the parameters of the five planets 
which have been found in this sample are given in Table [2j We also bring two independent 
simulation algorithms to bear on these data to derive extremely tight limits on planetary 
companions in the range 2d < P < 300d. The planetary periods and masses probed by these 



Jones et al. 


(2003 


) and 


Udrv et al. 


(2003) 



a lack of planets with orbital periods of 10-100 
days. Second, the large amount of high-precision data used here enables the detection of 
potentially terrestrial-mass planets. Core-accretion simulations which include Type I migra- 
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tion (llda &: Lid l2004j . l2008bl ) predict an abundance of planets with masses less than about 
10 M e orbiting within 1 AU, and a lack of planets with higher masses in that region (the 
"planet desert"). The large, extremely high-precision data set analysed here represents one 
of the first able to test these predictions. 

In Section 2, we describe the two methods used to derive detection limits from these 
data. Section 3 presents the results and discusses differences between the two methods. In 
Section 4, we further discuss the implications of our results with regards to the underlying 
distribution of exoplanets. 



2. Computational Methods 



There are two ways to appr oach the input data for a detection-limit c omputation: using 
the actual radial- velocity data ( Wittenmyer et al. 2006 ; Endl et ahlbood ). or creating simu 



lated data sets using the obs ervation dates and uncertainties of the real data (lO'Toole et al. 



2009al ; ICumming et al.lll999l ). We have employed independent methods using each of these 



approaches to determine the types of planets to which these data were sensitive. 



2.1. Simulated Data 



The first simulation technique, described fully in lO'Toole et al.l (j2009al ). seeks to re- 
cover simulated Keplerian signals (with appropriate noise terms arising from measurement 
uncertainties and stellar jitter) sampled at the epochs of the actual data and retaining the 
original uncertainties of the actual data. The si mulated velocity curves are analysed using 
the two-dimensional Lomb-Scargle periodogram (10 'Toole et al.l 120071 . l2009al ) which extends 
the standard Lomb-Scargle periodogram to search eccentricity as well as period space. 

The simulated observations for this study were constructed in the followin g way. We used 



a grid o f periods, eccentricities and planet masses as input in the same way as lO'Toole et al. 



(j2009a[ ). with periods in the range 2-320 days, eccentricities from 0.0 to 0.6 in steps of 0.1 
and masses of 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.5 and 1.0 Mj. Each set of parameters was 
simulated 100 times, leading to 72800 simulated velocity curves for each of the 24 stars. 

This approach does not use the standard F-test to determine the relia bility of a detec- 
tion. Instead, it uses four detection criteria as follows ( lO'Toole et al.ll2009al ). First, the RMS 
of the simulated observations must be greater than or equal to the RMS of the residuals of 
the best-fitting model to those data. That is, the scatter should decrease after the subtrac- 
tion of a Keplerian model. Second, the measured semi- amplitude must be greater than or 
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equal to twice the measured semi-amplitude uncertainty plus the RMS of the residuals of 
the best-fitting model. Third, the measured period must be greater than or equal to twice 
the measured period uncertainty. Finally, the value of x\ must be less than or equal to 
three. The la s t criter ion is only useful when stellar jit t er is n ot included, as discussed by 
O 'Toole et al.l ( ]2009al ); in this paper and lO' Toole et al.l ( l2009ch - where jitter is included - 
it is somewhat redundant. The other three criteria are a gauge of the reliability of a fit. 

With these criteria, the method does not rely on the assumptio ns of the F-test. Its 
chief strengths lie in its ability to recover significantly eccentric orbits ( jO'Toole et al.ll2009ai ) 
which would be undetectable using a standard Lomb-Scargle periodogram. In addition, 
because this method performs a full Keplerian fit to the simulated data, the distribution of 
orbital parameters derived fro m the fits provides va luable information on the biases inherent 
in fitting radial- velocity data (jO'Toole et al.ll2009al ). For the rest of this paper we will refer 
to this algorithm as Method 1. 



2.2. Real Data 



The second technique, described fully in IWittenmyer et al.l (120061 ) . adds the Keplerian 
signal of a simulated planet to the original data, and then at tempts to recover that signal 
using a Lomb-Scargle periodogram ( Lomb 1976[ Scargle 1982h . P revious work using similar 
methodology includes I Walker et al.l ( 119951 ) and lEndl et al.l (120021 ) . We will refer to this algo- 
rithm as Method 2. In this work, we required that the recovered period be within 5% of the 
injected period, and that its false-alarm probability (FAP) be less than 0.1%. As the null 
hypothesis of this method is that the dat a are purely nois e and contain no signals due to 
real pl anets, the known planets HD 4308b JlJdrv et al.ll2006h and HP 16 417b Jo 'Toole et al. 



2009bl ) must be treated differently. As noted in IWittenmyer et al.l ( 120091 ). the velocity signa- 
ture of an unknown planet can be altered by the presence of an additional (known) planet. 
H ence, we have determined the detection limits for HD 4308 and HD 16417 using the method 
of IWittenmyer et al.l (120091 ) wherein the simulated planet is added prior to fitting for and 
removing the known planet's signal. For the Alpha Centauri system (HD 128620/1), we 
removed a quadratic fit to the data, which provided a satisfactory approximation of the 
binary orbit over the short time period considered here. We also removed linear trends from 
the long-period binaries HD 10360 and HD 10361. Each star received 7 trials at the 99% 
recovery level, for trial signals with e =0.0, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6. In addition, 
simulations at 90%, 70%, 50%, and 10% recovery were performed at e = 0.0 for each star. 



This method is computationally faster than Method 1, since each trial only requires the 
computation of a simple periodogram, and there are considerably fewer such trials. Method 2 
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also has the advantage that the noise characteristics of the real data are preserved, whereas 
for simulated data, additional assumptions must be made about the noise term. However, 
it breaks down when there are fewer than N ~30 data points, as the reliability of the FAP 
calculation is strongly dependent on the number of data points. With limited data, even a 
very strong periodogram peak may fail to reach the 0.1% FAP level. This is less relevant for 
the curre nt work, as all 24 stars considered here have at least 31 observations (Tabled]). We 
note that ICumming et al.l (120081 ) in their Appendix give formulae for the FAP calculation in 
the regime of small N. For future work, the FAP calculations in Method 1 can be modified 
in this manner to better handle small N. 



3. Results 

3.1. Detectabilities 

For the purpose of this study, we define "detectability" as the fraction of simulated 
planets which were successfully recovered by the detection criteria employed by each method. 
To obtain a "big-picture" view of the types of planets detectable with our current data 
set, we have averaged the results over all stars and all eccentricities. Figure [T] shows the 
detectabilities from Method 1; the contours indicate the fraction of simulated planets that 
were recovered. The 99% detection limits obtained in this manner from Method 2 are 
overplotted as a dashed line. Results for each individual star are shown in Figures l2TfT3l 
The detectabilities obtained in this work are highly dependent on the characteristics of the 
data set; HD 20794 (Figure^]) has a large number of observations with little scatter (N = 78, 
rms=3.43 ms _1 ), whereas HD 28255a (Figure[6]) has few data points and considerable scatter 
(N = 37, rms=7.57 ms" 1 ). These results highlight the need to perform detection-limit 



computations on a star- by-star basis (jO'Toole et al.ll2009al ) rather than "whole of sample. 



Table [3] summarises the results from both methods for each star. We list the mean velocity 
amplitude K (averaged over all periods 2-300 days) which was recovered at the 99% and 
50% levels for each method. 



3.2. Comparing the Methods 

In general, the results from the two methods are in agreement, though for P < 10 days, 
Method 1 gives mass limits higher than Method 2 by a factor of 2-4 (cf. Fig. [TJ. The reason 
for this difference is that Method 1 incorporates the added criterion that the semi-amplitude 
K of the recovered Keplerian orbit must satisfy K/ctk > 2. Removing this criterion results 
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in an extremely high proportion of false positives (incorrect detections). Method 2, with a 
periodogram recovery, has no such constraint on the amplitude, only on the period. Both 
are more conservative than a human would be in analysing an individual observed data set 
by hand. On the other hand, they have a quantitative selection process which means that, 
given their detection criteria and a set of data, one can use these algorithms to draw robust 
conclusions about planetary frequency and orbital parameter distribution. 

Both Methods 1 and 2 address the question of "What planets could have been detected 
by these data?" In any such detection-limit determination, the aim is to automate and 
perfor m thousands of times what a human investigator would reasonably do to detect a 
signal (10 'Toole et al.ll2009af ). This process necessarily involves some sacrifices. For exam- 
ple, Method 2 attempts to detect an injected signal using the standard Lomb-Scargle peri- 
odogram, an extremely fast computation. However, the Lomb-Scargle periodogram is most 
efficient at detecting sinusoidal signals; when Kepleria n orbits become moderately eccentric. 



their shapes can become signifi cantly non-sinusoidal (ICumming et al. 



2008 



O'Toole et al. 



2009al ; IWittenmyer et al.ll2009l ). To work around this issue, IWittenmyer et al. 



re- 
stricted the eccentricities of the injected planets to e < 0.6, a range which includes about 
90% of known exoplanets. The effects of eccentricity remain evident, however, in the "blind 
spots" where Method 2 indicates that planets at certain periods are completely undetectable 
(Figure [T5|) . For each star, there are certain trial periods where the data sampling conspires 
with an injected eccentric signal, such that an ordinary Lomb-Scargle periodogram does not 
exhibit a significant peak at that period. Method 2 would then reject the trial signal, even at 
a large amplitude (K max = 100m s _1 ), and these situations result in allegedly undetectable 
trial signals. Method 1, however, using the 2DKLS periodogram which is much more effective 
at detecting eccentric signals, would not fall victim to this pathology. 

The significance of a periodogram peak obtain ed by Method 2 is computed using the 
analytic formula given in iHorne fc Baliunasl (119861 ). However, the false-alarm probability 
(FAP) obtained by this formula tends to b e smaller (i.e. more s ignificant) than that obtained 
by a bootstrap randomisation procedure (jKiirster et al.l 119971 ). This is due to the fact that 
noise in radial- velocity data tends to be non-Gaussian. Thus, a human would visually inspect 
the periodogram, then perform a bootstrap FAP calculation to obtain an accurate estimate 
of the believability of any peaks. The bootstrap FAP method, however, requires hours to 
calculate, rather than the milliseconds it takes for the analytic method, and can therefore 
not be usefully employed in the analysis of tens of thousands of simulations. 

Another necessary concession is the requirement that the recovered period be within 5% 
of the injected period. The purpose of this constraint is to exclude significant periodogram 
peaks at aliases of the true period. Spurious periods were also often produced by the lunar 
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cycle (~30 days), and the duration of the densely s ampled "Rocky Planet Search" (48 d ays). 
All previously published work using this method ( 1 Wittenmyer et al.l 120061 . 120071 120091 ) has 
used this 5% criterion. This has the effect of preventing incorrect detections (false positives). 
Method 1 does not have an explicit criterion of this sort, but the four criteria outlined in 
§ 2.1 serve to keep the rate of false posi tives below ~l-2%; a detailed discussion of false 
positives is given in lO'Toole et al.l ( 12009al ) . 

To summarise, Method 1 seems to be most appropriate for significantly eccentric orbits 
because the simulated planets are fitted with a full Keplerian model. Since a value of stellar 
velocity noise (jitter) is added to the simulated data in Method 1, that method is more 
effective when the radial-velocity noise is well-understood. If the stellar jitter is not well- 
understood, it is more appropriate to use Method 2 since the actual velocity data (which 
include all noise sources) are used. 



4. Discussion 



4.1. M„ >100 M fl 



Jones et al.l (120031 ) and lUdry et al.l (120031 ) noted an emerging "period valley" between 
10 and 100 days; t he ^ 250 additional planets discovered since then have continued this 



trend. lUdry et al.l (120031 ) also showed that this feature is mainly attributable to the lack 



of massive (m sin i > 2Mj up ) planets in that period range. Our results support this claim, 
since the 24-star sample contains three low-mass planets with periods between 10-100 days 
(HD 4308b, HD 16417b, HD 115617c), and the data can exclude planets with m sin % > 
0.5Mj up at the 99% level for all of this region (Figure H}. The deficit of higher-mass planets 
with P < 300 days is readily apparent, suggesting that the period valley is no t a selection 
effect; this is consistent with the conclusions reached by ICumming et al.l (120081 ) . 



4.2. 



M p <100 M a 



At lower masses (10-100 M©) however, the situation is far less clear. This is an extremely 
interesting mass range, as core-accretion theory, and numerical simulations based on it, 



(Ida & Lin 


2004. 


2008a: 


Liu et al. 


2009) 



10Mg 

showed a pile-up of planets at 0.2-0.3 AU (P ~30-50 days). This is a mass range, then, 
which should show detectable features in the structure of the period distribution. However, 
this is also a mass range at which serious selection effects impact on the detectability of low- 
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mass planets, a nd most especially in the critical 15-60d orbital period range. This is why 



Liu et al.l ( 120091 ) noted that the predicted pile- up had not yet been observed. Figure [161 shows 
the period distribution of planets with M p <100 M e . A steep fall-off in the distribution is 
evident at P > 5 days. This can be understood as arising from two selection effects. First 
and most obviously, longer-period planets have weaker radial- velocity signals. Second, planet 
search observing runs are almost exclusively scheduled on large telescopes during bright time, 
resulting in a large peaks in the window function at 28d. The only means to mitigate this 
sampling problem is to observe through dark lunations. This is a strategy that only the 
AAPS is currently employing in its "Rocky Planet Search" observing campaigns and which 
(as demonstrated by the discovery of 61 Vir c, P = 38d, K = 3.6 ms" 1 ) has proven successful 
in finding low-mass planets in this period range. 

When observing strategies are employed which mitigate selection effects in this period 
regime, we do indeed find low-mass planets (e.g. Vogt et al. 2010, O'Toole et al. 2009b). These 
discoveries indicate that the observed fall-off in the distribution of such planets (Figure [T6|) 
may - contrary to what is seen at higher masses - indeed merely be a result of observational 
biases. If we look at the results of our simulations, they indicate that 10% of planets with 
K ~2.7 ms -1 could have been detected in our sample of 24 stars. One such planet with an 
even smaller radial-velocity amplitude has atually been discovered in our sample (61 Vir b, 
K = 2.2 ms -1 ), suggesting that at least a further nine more such planets could have been 
missed, and will remain below the detection threshold until additional data are obtained. 

In addition, the solar-type dwarf stars typically observed by precision radial-velocity 
programs have typical rotation periods in the 20-50 day range. Hence, emergent periodicities 
are not usually believed until sufficient data are in hand to be confident that the periodicity is 
planetary in origin rather than due to rotational modulation by starspots. This selection bias 
suggests that low-mass planets with P ~30-50 days may be present in extant radial-velocity 
data sets, but as yet not believed and unpublished. 



4.3. The underlying planetary mass function 



1980; 



he importance of plane tary migration has been appreciated for some time (Goldr eich fc Tremaine 



:e oi p 

Lin fc Papaloizoulll986l ). The modelling work of lTrilling et al.l ( 120021 ) and lArmitage et al 



(120021 ) has provided evidence that observational data for exoplanets are consistent with a 
model where planets are formed at around 5 AU and then undergo migration. This leads to 
a distribution of planetary semimajor axis a that depends on log(a). Although this scenario 
provides a good match to observations (e.g. Armitage 2007), migration rates at s hort orbital 



periods are expected to be different for different mass planets. The simulations of iRice et al. 
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( 120081 ) indicate three regimes: (1) low-mass planets (<0.5Mj up ) migrate inwards and stall 
at approximately half the rotation period of the star, (2) high-mass planets (10M Jup ) should 
collide with the star, and (3) the intermediate- mass (lM, Tlip ) plane t population should have 



an inner edge located at twice the Roche-limit (IRasio & Fordlll996l ). For the purposes of our 
models, we implement a distribution of planetary semimajor axis that depends on log(a), 
referred to here as "logarithmic migration." 

Although our target star sample and number of detections is small, we can use the fact 
that it is well-characterised in an attempt to understand the featu res of the large r inhomoge- 



nous data set of radial velocity exoplanets. Following the work of lOToole et al.l f2009c) and 
using our computed detectabilities and a simplistic logarithmic migration (jArmitagel 2007 ). 



Figure [TTJ shows the expected number of planets for different underlying mass functions as 
a function of semi-major axis for the case of equal numbers of exoplanets in each bin. The 
bold asterisks represent the planets detected from this sample (Table [2]). We find that the 
total expected number of detections for different assumed mass functions observed with our 
sample varies by a factor of about two for the mass functions shown, with all giving a smooth 
fall-off in the number of detections toward longer periods. Thus, although the peak of the 
histogram serves to approximately characterise our detection of two exoplanets in the 10-20 
day period bin and is consistent with those below the threshold, the clear rise at longer peri- 
ods in the observed distribution in Figure [17] is not reproduced. However, it must be noted 
that the sample of all radial velocity exoplanets shown in Figure [17] is highly inhomogeneous 
and represents an observed sample of a higher mass range. 

Although Figure [T7J indicates that the distribution in period space is relatively insensitive 
to different mass functions, we find greater sensitiv ity when detectability is plotted as a 



function of mass. Following lO'Toole et al.l (l2009d ). we show in Figure [18] the expected 



detections for our sample as a function of mass for a variety of assumed mass functions. 
We find that the peak in the expected number of observed exoplanets moves for different 
mass functions, a—— 1.7 leads to a peak in the number of detected planets around 0.02 Mj up 
(6.3 M ffi ), the a=— 1.2 expected detections peak at around 0.05 Mj up (16 M ffi ) while the flatter 
a=— 0.7 mass function produces an expected detection peak at around 0.10 M Jup (32M ffi ). 

Given that we have found two exoplanet signals in the 16-25 M e minimum mass range 
(corresponding to the 0.05 Mj up or 16 M e bin in Figure [18] we may consider these as providing 
rough limits on the underlying mass function as determined by how changes in the mass 
function cause the peak of detections to move toward and away from the 0.05 M Jup (16 M e ) 
bin. We find the peak moves from the 0.02 to the 0.05 Mj up bin for a values shallower than 
a=— 1.6 and from 0.05 Mj up to 0.1 Mj up at a—— 1.1. It is noteworthy that the peak mass 
for different mass functions is sensitive to migration. If no migration is employed then the 
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mass function peak moves away from the 0.05 Mj up (16 M e ) bin when 0.6 < a < 1.2. Based 
on Poisson-counting statistics of two detections our s imulations indicate a normalisation of 
0.17±0.12 (corresponding to the F value presented in lOToole et al.l (12009a )). Our estimate 
of the unde rlying mass fu n ction shown in Figure [TS] is consistent with the first attempt 
presented in lO 'Toole et al.l ( ]2009d ). That is, the additional planet detections in this sample 
and the comprehensive detection limits presented here serve to strengthen the case for an 
underlying planet mass function of the form dN/dM oc M~ L0 . Given that we robustly do 
not detect more massive planets in this sample, flatter mass functions are not favoured. 



5. Conclusions 

We have shown, using two independent simulation algorithms, that recent data from the 
AAPS are capable of placing meaningful limits on the population of potentially terrestrial, 
sub-Neptune-mass planets orbiting nearby stars. As we are able to have detected all planets 
with at least a Saturn mass in the period range P < 300 days, we conclude that the observed 
lack of massive planets (M p >100 M ffi ) with periods 10-100 days is real and not due to 
selection effects. However, we conclude that the shortage of lower- mass planets (M p <100 
M e ) predicted by core-accretion theory in this period range may arise from observational 
selection biases against detecting such planets. 
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Table 1. Summary of Radial- Velocity Data 



Star 


N 


RMS 
(m s _1 ) 


< cr > 
(ms- 1 ) 


HD 1581 


58 


2.81 


0.95 


HD 4308 


77 


2.92 a 


1.09 


HD 10360 


31 


3.77 b 


0.90 


HD 10361 


31 


3.20 b 


0.84 


HD 10700 


95 


2.65 


0.86 


HD 16417 


55 


2.66 a 


0.86 


HD 20794 


78 


3.43 


0.81 


HD 23249 


64 


3.16 


0.56 


HD 26965 


55 


4.80 


0.76 


HD 28255A 


37 


7.57 


1.18 


HD 43834 


78 


4.27 


0.80 


HD 53705 


80 


3.09 


1.21 


HD 72673 


50 


2.66 


1.06 


HD 73524 


44 


4.22 


1.28 


HD 84117 


83 


5.00 


1.33 


HD 100623 


62 


4.08 


1.00 


HD 102365 


97 


2.49 


0.89 


HD 114613 


106 


4.25 


0.78 


HD 115617 


94 


4.29 


0.74 


HD 122862 


49 


3.81 


1.27 


HD 128620 


33 


1.94 c 


0.40 


HD 128621 


42 


2.76 c 


0.56 


HD 136352 


78 


4.34 


1.07 


HD 146233 


48 


3.91 


1.02 



a After removal of known planet's orbit. 
b Residuals of linear fit. 
c Rcsiduals of quadratic fit. 



Table 2. Planets From This Sample 



Planet 


Period 


To 


e 




K 


M sin i 


(4 


Reference 




(days) 


(JD-2400000) 




(degrees) 


(ms _1 ) 


(M ffi ) 


(AU) 




HD 16417 b 


17.24±0.01 


50099.7±3.3 


0.20±0.09 


77±26 


5.0±0.4 


22.1±2.0 


0.14±0.01 


O'Toole et al. (2009b) 


HD 4308 b 


15.609±0.007 


50108.5±1.9 


0.27±0.12 


210±21 


3.6±0.3 


13.0±1.4 


0.118±0.009 


O'Toole et al. (2009c) 


HD 115617 b 


4.2150±0.0006 


53369. 166± (fixed) 


0.12±0.11 


105±54 


2.12±0.23 


5.1±0.5 


0.050201±0.000005 


Voet et al. (2010) 


HD 115617 c 


38.021±0.034 


53369. 166± (fixed) 


0.14±0.06 


341±38 


3.62±0.23 


18.2±1.1 


0.2175±0.0001 


Voet et al. (2010) 


HD 115617 d 


123.01±0.55 


53369. 166± (fixed) 


0.35±0.09 


314±20 


3.25±0.39 


22.9±2.6 


0.476±0.001 


Voet et al. (2010) 
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Table 3. Summary of Detection Limits at e=0.0 







Mean K velocity recovered (ms 


V ) 




99% recovery 


99% recovery 


50% recovery 


50% recovery 


O LeLI 




Method 2 


Method 1 


A/Tot Vinrl 9 


HD 1581 


12.6±1.1 


4.8±2.2 


6.0±0.5 


2.8±0.6 


HD 4308 


5.9±0.1 


4.7±2.0 


3.2±0.1 


2.7±0.8 


HD 10360 


11.9±0.3 


9.6±1.8 


3.4±0.1 


7.0±3.1 


HD 10361 


14.2±0.3 


7.9±2.0 


3.4±0.1 


5.6±2.5 


HD 10700 


5.0±0.2 


3.5±0.9 


2.8±0.1 


2.4±0.4 


HD 16417 


5.8±0.3 


5.9±2.5 


3.1±0.2 


3.4±0.9 


HD 20794 


4.8±0.1 


4.3±1.0 


2.4±0.1 


2.7±0.5 


HD 23249 


12.7±3.0 


4.6±1.9 


6.2±1.4 


2.8±0.7 


HD 26965 


5.3±0.6 


8.3±3.7 


2.5±0.3 


6.2±1.8 


HD 28255A 


34.2±10.2 


14.1±3.2 


8.8±2.4 


11.4±5.5 


HD 43834 


4.6±0.4 


6.3±2.1 


2.5±0.2 


4.0±0.8 


HD 53705 


5.7±0.3 


5.0±2.6 


3.1±0.2 


3.0±0.9 


HD 72673 


14.5±0.5 


5.4±2.8 


3.3±0.1 


4.7±4.3 


HD 73524 


14.4±1.2 


8.6±3.8 


4.1±0.3 


6.8±4.0 


HD 84117 


7.9±0.1 


7.9±3.4 


4.1±0.1 


5.1±1.5 


HD 100623 


6.5±0.8 


7.2±3.5 


2.9±0.4 


4.8±1.5 


HD 102365 


5.6±0.1 


3.9±1.6 


2.9±0.1 


2.5±0.6 


HD 114613 


8.3±0.1 


6.7±2.9 


4.3±0.1 


4.5±1.0 


HD 115617 


4.0±0.2 


7.3±3.2 


2.4±0.1 


4.8±1.4 


HD 122862 


8.9±1.1 


6.4±2.0 


4.0±0.5 


4.8±1.7 


HD 128620 


19.4±2.3 


4.8±2.3 


3.1±0.3 


4.1±3.6 


HD 128621 


16.7±0.7 


8.6±8.1 


3.1±0.1 


6.0±4.5 


HD 136352 


6.0±0.4 


7.4±2.8 


3.1±0.2 


5.0±1.4 


HD 146233 


6.2±0.1 


6.9±2.6 


3.2±0.1 


4.2±1.2 
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10 100 

p (d) 

Fig. I. — Detectabilities obtained from Method 1 (solid contours) averaged over all 24 stars, 
for planets with eccentricities from 0.0 to 0.6. Contours indicate the fraction of injected 
planets that were recovered. Results from Method 2 (99% recovery) for the same range of 
eccentricities are shown as a dotted line. The five announced planets in this sample are 
plotted as large asterisks. 
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Fig. 2. — Detectabilities obtained from Method 1 (solid contours), for HD 1581 (left panel) 
and HD 4308 (right panel), for planets with eccentricities from 0.0 to 0.6. Contours indicate 
the fraction of injected planets that were recovered. Results from Method 2 (99% recovery) 
for the same range of eccentricities are overp lotted as a dotted line. 
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;. 5.— Same as Figure El but for HD 20794 (left) and HD 23249 (right). 
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;. 6.— Same as Figured but for HD 26965 (left) and HD 28255A (right). 
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7.— Same as Figure but for HD 43834 (left) and HD 53705 (right). 
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Fig. 8.— Same as Figure [2 but for HD 72673 (left) and HD 73524 (right). 
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Fig. 9.— Same as Figured but for HD 84117 (left) and HD 100623 (right). 
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Fig. 10.— Same as Figure El but for HD 102365 (left) and HD 114613 (right). 
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;. 11 — Same as Figured but for HD 115617 (left) and HD 122862 (right). 





P (days) 



P (days; 



;. 12.— Same as Figure but for HD 128620 (left) and HD 128621 (right). 
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;. 13.— Same as Figure but for HD 136352 (left) and HD 146233 (right). 
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Fig. 14. — Detection lim its for HP 115617 (61 Vir), which has recently been found to host 
three low- mass planets (jVogt et al.l 120101 ) . These planets were not accounted for in the 
detection-limit simulations presented here. The lines represent recovery rates of 99% (solid), 
50% (dashed), and 10% (dotted). The planets are shown here as filled circles. The discovery 
of these planets was only possible with the inclusion of data from the Keck telescope and 
from a second 47-night Rocky Planet Search campaign in 2009 July/August. In this work, 
we consider data from the AAT only, from 2005 June-2009 June. 
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Fig. 15. — Detection limits for HD 1581 (Method 2) for test signals with e = 0.0 (solid 
line) and e = 0.6 (dashed line). At higher eccentricities, there arise some periods for which 
the test signal is never recovered at sufficient significance by the traditional Lomb-Scargle 
periodogram. 
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1 10 100 1000 

Period (days) 

Fig. 16. — Observed orbital period distribution of radial- velocity planets with m sin % <100 
M e (N=50). The fall-off at P > 5 days arises from the twin selection effects of smaller 
velocity amplitudes and poorer temporal sampling, especially in the period range of 30-50 
days. 
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Fig. 17. — Solid-line histogram: Orbital period distribution of radial- velocity exoplanets 
with m sin % > 0.5Mj up (N=256); both simulation methods discussed in this work show that 
all of the planets this massive would have been detected in our sample, indicating that the 
lack of giant planets with P < 100 days does not arise from selection biases. Overplotted 
as dashed and dotted histograms are orbital period distributions that would be expected 
from our simulations, for a variety of assumed power-law mass functions. Asterisks show 
the periods of announced exoplanets from our sample that are redetected by our automated 
selection criteria (§ 2), while arrows show the periods of the 3 planets detected orbiting 61 Vir 
(using a combination of AAPS and Keck data), which are not redetected by our automated 
criteria with AAPS data alone. 
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Fig. 18. — Expected number of exoplanets from our sample as a function of mass, for 
different mass functions including logarithmic migration. Asterisks and arrows show the 
mass of announced exoplanets from our sample, as for Figure fT71 The heavy filled circle with 
uncertainties shows the number density represented by our automated system detections in 
the 0.05Mj up bin, while arrows show the periods of the 3 planets detected orbiting 61Vir 
(using a combination of AAPS and Keck data), which are not redetected by our automated 
criteria with AAPS data alone. 



